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Abstract 

A cosmology with the growing cosmological term is considered. If there is no 
exchange of energy between vacuum and matter components, the requirement of 
general covariance implies the time dependence of the gravitational constant G. 
Irrespectively of the exact functional form of the cosmological term growth, the 
universe ends in a de Sitter regime with a constant asymptotic A, but vanishing 
G. Although there is no divergence of the scale factor in finite time, such as in the 
"Big Rip" scenario, gravitationally bound systems eventually become unbound. 
In the case of systems bound by non-gravitational forces, there is no unbounding 
effect, as the asymptotic A is insufficiently large to disturb these systems. 



Cosmological observations of increasing quantity, quality and diversity have 
established a new picture of the universe, its composition and dynamics pp. Mea- 
surements of luminosity-redshift relations for the supernovae of type la (SNIa) [2] 
and the temperature anisotropies of the cosmic microwave background radiation 
(CMBR) 3 , in agreement with other cosmological observations, have determined 
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that the universe is presently in the phase of accelerated expansion, which is at- 
tributed to a new component of the universe named dark energy. The origin of 
dark energy (which together with dark matter constitutes the "dark sector" of the 
universe) has not yet been clarified. Pretenders to the title of dark energy are, 
however, numerous. The first in a long line of candidates is certainly the cosmo- 
logical constant (CC), also known as cosmological term (UG3E1- This conceptually 
simple choice, which naturally fits into the formalism of general relativity (GR), is, 
however, burdened with problems which, to be solved, require a highly unnatural 
amount of fine-tuning. Namely, quantum field theory contributions to the vacuum 
energy (which is equivalent to the CC) and the observed value of the CC differ 
by many orders of magnitude. Furthermore, there is a coincidence problem of 
understanding why the CC energy density and the non-relativistic matter energy 
density are comparable at the present epoch. These problems were a strong incen- 
tive towards the development of dynamical dark energy models. One of the most 
studied classes of such models are quintessence models [Zj. Quintessential models 
describe dark energy in terms of the scalar field slowly rolling in a potential. An- 
other appealing proposal is the Chaplygin gas 0. The most interesting feature of 
this model is the possibility of unified description of the "dark sector" , i.e. the uni- 
fication of dark matter and dark energy. Another alternative to quintessence with 
the similar properties of effective unification of dark sector is Cardassian cosmol- 
ogy 0- Other interesting models include the non-perturbative effects of vacuum 
energy jlUj and the fluctuating CC jllj . The consideration of the CC as a dynam- 
ical quantity in the framework of the renormalization group equation |12[ I13| 11 4 j 
is another promising approach. 

A large majority of dark energy models describes dark energy in terms of the 
equation of state (EOS) 

Pd = wp d , (1) 

where w is the parameter of the EOS, while pd and pd denote the pressure and 
the energy density of dark energy, respectively. The value w = —1 is characteristic 
of the cosmological constant, while the dynamical models of dark energy generally 
have uu > — 1. A recent examination of the dark energy EOS, based on the data 
from CMBR, SNIa, large scale structure (LSS) and Hubble parameter measure- 
ments from the Hubble Space Telescope (HST), assuming the redshift independent 
parameter w, restrict w to be in the interval —1.38 < w < —0.82 at the 95% 
confidence level |15| . This result raises a question about the possibility of dark 
energy models with the supernegative EOS, i.e. with w < —1. This new sort of 
dark energy, first analyzed in which was followed by numerous analyses |17j . 
has soon deserved a term of its own, phantom energy. The most interesting feature 
of phantom energy models is the possibility of a divergence of the scale factor of 
the universe a in finite time. Such a behaviour of the scale factor has a dramatic 
effect on all bound systems. Namely, the bound systems become unbound at some 
finite time interval before the onset of the divergence in a. This type of the fate of 
the universe is known as a "Big Rip" scenario |18j . 

In this letter we consider the case of the growing cosmological term A and its 
implications for the asymptotic expansion of the universe and the destiny of the 
bound systems. Let us start with the specification of the main characteristics of 
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the components of the universe. We assume that we have two components of the 
universe: non-relativistic matter and the variable cosmological term. The inclusion 
of radiation or other components is straightforward, but without any substantial 
influence on the future evolution of the universe. Non-relativistic matter has the 
EOS 

Pm = 1pm , (2) 

where 7 > is the parameter of the EOS and p m and p m represent the pressure 
and the energy density of non-relativistic matter, respectively. The conservation 
of the energy-momentum tensor of non-relativistic matter, Tj^ v = 0, leads to the 
standard relation of the evolution of p m with the scale factor 

/ a \ - 3 ( 1+ t) 

Pm = Pm,0 — • (3) 

The equation of state of the cosmological term is 

PA = ~PA , (4) 

where p\ is the pressure and p\ is the energy density of the cosmological term. The 
requirement of the conservation of the energy-momentum tensor of cosmological 
term, T^" = pa9^ u , leads to a constant cosmological term energy density p\. If we, 
however, consider a variable (time-dependent) p\, T^ u can no longer be conserved. 
Models with variable A = 8nGp\ and G were extensively studied in The 
consideration of phantom energy models in an analogous framework leads to an 
interesting general finding on the future evolution of the universe |2U| . In the 
setting described above, the general covariance of the Einstein equation 

Gy, u = -SttGT^ , (5) 

where G^ is the Einstein tensor and T^ v = + Tfc u , can be maintained if 
the gravitational constant G acquires space-time dependence. This fact can be 
interpreted as a modification of the dynamics of General Relativity. We effectively 
describe this additional dynamics by promoting G into a function of space-time 
coordinates. We further assume that G depends on time only, i.e. G = G(t) 1 . 
The condition {G{t)T^ u )- v = can be expressed as 

d(G( Pm + p A )a 3 ) = -G(p m + p A )da 3 , (6) 

which gives the law of evolution of G(t) 

G( Pm + pa) + Gp A = . (7) 

Apart from the procedure displayed above, there exist other ways how the vari- 
ability of the cosmological term can be incorporated into the formalism of General 
Relativity. One possibility is to allow the interchange of energy and momentum 
between the cosmological term (vacuum) and matter and radiation components 

interesting models with variable G are also obtained in the formalism of Quantum field theory in 
the curved space-time |21|. 
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|22j . In this model, the energy-momentum tensors of separate components (non- 
relativistic matter or radiation and vacuum energy) are not conserved, but the 
total energy-momentum tensor is. Furthermore, G is time-independent in this set- 
ting. It is interesting to see that both the model of our paper and the model given 
in |22| follow from equation © under different assumptions. The dynamics of our 
model © is obtained assuming that the vaccum energy is time dependent and that 
the tensor of energy-momentum of the matter component is separately conserved. 
The dynamics of [22], expressed by their equation (5a), is recovered under the as- 
sumption that G is time-independent and that there exists an exchange of energy 
between the vacuum and the matter components. 

The set of equations governing the evolution of the universe is completed by 
the Friedmann equations for the scale factor 

d\ 2 k 87r . . 

- + — = -Z-G(p m + P A , 8 
a J a z 6 

- = — 7rG{p m + p\ + 3p m + 3p\) . (9) 
a 3 

Given the set of equations ©, ©, © and ©, we can investigate the future 
evolution of the universe for a general growing p\. The law of evolution of p m © 
clearly shows that at sufficiently distant future times we have p A p m . In this 
limit, equations © and © become 

Gp A + Gp A = 0, (10) 

'a\ 2 k _ 8vr 
, a / a 2 3 

From IjlOJI it follows that at distant times we have Gp A = const. Therefore, in 



+ ^2=— GPA- (11) 



equation (jll|) we can disregard the k/a 2 term at sufficiently distant times. We 
finally obtain the equations governing the evolution of the universe 

d(G PA ) 



dt 

a\ 2 8vr 



0, (12) 



, ~rrGp\ . (13) 
a) o 

The solution of the equation displayed above is of the form 

a ~ e V^F* ; (14) 

where Aoo = 87r(G ? pA)t-*oo- in the distant future the universe enters the de Sitter 
regime of expansion. This fact does not depend on the concrete form of p\, as long 
as it is a growing function of time (explicitly or implicitly via some other quantity 
depending on time, such as the scale factor a). Furthermore, the growth of p A 
implies that G is a decreasing function at distant times. The unbounded growth 
of p A results in vanishingly small values of G at distant times. 

Let us further examine the complete evolution of the universe for a class of 
growing p A models. We assume the following form for the variable cosmological 
term energy density: 

/ a \ - 3 ( 1+r ?) 

PA = PA,0 — , (15) 

\ao/ 
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where rj < — 1. We consider the flat universe case, k = 0, which implies that 
the total energy density po equals the critical energy density p Ci o, and introduce a 
dimensionless parameter Q A = PA,o/po- The subscript or superscript denotes the 
present time throughout the paper. What remains is solving equations (0) and (jHJ) 
given the laws of evolution of the energy densities © and (|T5)) . The quantity G 
evolves according to the law 



l + V 

v — 1 



G = G [n° A {-) 3(v y) + i-n° A ) tin) 



while the scale factor is implicitly given by the expression 

/a/aa , *. 11+2. 

X 2(i+37) (nlx- 3 ^) + i _ j 2 "" 7 dx . (17) 

The evolution of the scale factor of the universe a with time is shown in figure 
n The dependences of a on time for different values of the parameter 77 differ 
from each other the most in the distant past and the distant future. Graphs of the 
evolution of the scale factor in the distant future also reveal the beginning of the 
exponential expansion, i.e. the onset of the de Sitter regime. The time dependence 
of G(t), depicted in figure is more sensitive to the value of the parameter 77. 
Figure EJ displays some general features of the dynamics of G(t). In the early 
universe (for small values of the scale factor) the change of G(t) is very slow. The 
rate of change is the greatest around the present epoch, while at large times G(t) 
tends towards 0. The pronounced dynamics of G(t) for more negative values of 
the parameter rj would probably be the best testing ground for the proposed class 
of models, i.e. it would provide the most stringent constraint on the growth of the 
cosmological term energy density with the scale factor, described by the parameter 
77. The time evolution of A is displayed in figure El In the early universe the value 
of A approaches for 77 < — 1. For large times, the function A(i) tends to its 
asymptotic value Aqo. For more negative values of 77, the asymptotic value A^ 
increases and it is approached slowlier. The dependence of the asymptotic value 
Aqo on the parameter 77 is shown in figure |IJ For the chosen range of the parameter 
77, approximately consistent with the constraint of reference |15j . the quantity Aqq 
grows quite modestly, definitely remaining of the same order of magnitude. 

The evolutions of A(t) and G{t) depicted in figures and 03 respectively, show 
that the universe ends up in a state with finite Aoo and vanishing G^. An impor- 
tant question for such a cosmology is the destiny of bound systems. First, we focus 
on gravitationally bound systems. Let us consider a spherically symmetric system 
with a mass M at the centre, in the cosmology with A^ ~ const and G^ ~ 0. 
We can treat this system as approximately static. The metrics of such a system is 
then given by [23] 

ds 2 = e v dt 2 - e x dr 2 - r 2 dQ 2 , (18) 



where 



A 2GooM l.o 

e~ X = 1 -A^r 2 . 19 

r 6 



In the Newtonian limit, the gravitational potential becomes 

G00M 1 2 



(r) = — -rAooT^. (20) 

r b 
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As Gqo is vanishingly small and Aoo is finite, presently gravitationally bound sys- 
tems cannot remain bound. The unbounding happens at some instance before a 
state with A»o w const and G^ « is achieved. Clusters of galaxies, galaxies, 
stars and stellar systems will fall apart. Planets will lose their atmospheres. The 
gravitational interaction will become fully dominated by the asymptotic value of 
the cosmological term. 

Next, we turn to systems bound by forces other than gravitational, e.g. strong 
or electromagnetic. In these systems, G^ is clearly of no relevance. The depen- 
dence of Aqo on the parameter n is shown in figure For values of r/ not much 
smaller than —1, Aqo remains of the same order of magnitude as Ao- As Ao obvi- 
ously has no unbounding effect on non-gravitationally bound systems, and given 
the difference of many orders of magnitude between typical interaction scales of 
gravitational and non-gravitational forces, neither A^ will be able to unbound 
non-gravitationally bound systems. The break-up of atoms and nucleons induced 
by the variable cosmological term will not happen. 

Finally, let us discuss constraints on the model parameters. Cosmological ob- 
servations impose rather strong constraints on the time dependence of the gravita- 
tional constant G. The observational bounds on the variation of G differ in their 
origin and refer to different epochs of the universe expansion |24| I25| . The observa- 
tions of spin-down rate of pulsars yield constraints \G/G\ < (2.2 — 5.5) x 10~ n yr -1 
for the pulsar PSR B0655+64 [II], and \G/G\ < (1.4 - 3.2) x lO^V -1 for the 
pulsar PSR J2019+2425 |25l I27j . where the range in constraints comes from the 
uncertainties in the neutron star equation of state. The effects of the variation of G 
on the orbital period of the neutron star- white dwarf binary system PSR B 1855+09 
give the constraint G/G = (-1.3 ± 2.7) x Kr n yr _1 [23 HH]- Tne consideration 
of effects of the variation of G on the Chandrasekhar mass in double-neutron- 
star binaries yields a constraint of G/G = (—0.6 ± 4.2) x 10 _12 yr _1 The 
observations of SNIa give G/G < 10 _11 yr _1 at redshifts z ~ 0.5 [301) while the 
Big Bang Nucleosynthesis constrains the gravitational constant during BBN to be 

Gbbn/Go = 1-01 ~^~q'iq a t the 68.3 % confidence level One of the strongest 

constraints on the variation of G comes from the Lunar Laser Ranging (LLR) and 
amounts to G/G = (0.46 ± 1.0) x 10~ 12 yr _1 32 J. Figure El shows that constraints 
on G/G from the distant past are more easily satisfied than those from the recent 
past since the variation of G is stronger in the recent than in the distant past. 
Namely, the rate of change G grows with time in absolute value from the distant 
past to the present epoch, while the value od G decreases with time. For given 
n, of all past epochs, the quantity G/G has the strongest variation for the present 
epoch. On the other hand, the strongest observational limit on G/G comes from 
LLR 32 at the present epoch (with z = 0). Therefore, the comparison of theo- 
retical expressions and observational bounds at the present epoch can produce the 
most stringent constraints on the parameter n. From equations Q and (|15|) the 
theoretical expression for the rate G/G at z = becomes 



Using the value Hq = 72 km/s/Mpc for the Hubble parameter, Q° = 0.7 and the 




(21) 
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values for G/G from jH2j, we obtain the following constraints on the parameter r\: 
i] > —1.0035 at the la level, rj > —1.01 at the 2a level and rj > —1.016 at the 3a 
level. Other observational constraints yield weaker limits on rj. The constraints 
from BBN [31], e.g., using (fTHjl yield 77 > —1.1881 at the 68.3 % confidence level, 
which is a much weaker limit. All the figures are plotted using a broader range 
of parameter r/ than the range constrained by LLR to better illustrate the form 
of dependencies of functions a, G and A on time and rj. The obtained constraints 
result in strong limits on the parameter rj which must be close to -1. However, as 
long as the condition r/ < —1 is satisfied, no matter how close 77 is to -1, the main 
conclusions of this paper remain unaltered. In another words, even for very slow 
growth of pa (which satisfies all the conditions on the variation of G) the partial 
rip scenario is valid, i.e., in the distant future the gravitationally bound systems 
become unbound, while the nongravitationally bound systems remain bound. 

Considerations presented above show that a universe with the growing cosmo- 
logical term and the time-dependent G has very interesting properties. The energy 
density of dark energy component grows with the scale factor (like in phantom en- 
ergy models), but the parameter of the EOS remains -1 (unlike in phantom energy 
models where the same parameter is less than -1). The universe ends up in a de 
Sitter regime with constant Aoo and vanishing Goo. The expansion continues in- 
finitely, without any divergence of the scale factor in finite time. However, despite 
the absence of the abrupt ending of the universe in a "Big Rip" event, accompanied 
by the unbounded growth of dark energy density, all bound systems do not remain 
bound. Gravitationally bound systems become unbound in the distant future in an 
interplay of constant repulsion of and vanishing attraction mediated by G^. 
Systems bound by non-gravitational forces, however, do not share the destiny of 
their gravitationally bound counterparts. Numerical calculations show that, for 
a reasonable range of the parameter 77, the value A^ is not sufficiently large to 
disturb systems bound by, e.g. strong or electromagnetic forces. It is important 
to emphasize once more that all these findings are not dependent on the size of 
77, as long as it remains less than -1. For the case of r] close to -1, all the effects 
described above happen, only the variations of G(t) and A(t) are much milder and 
the onset of the saturation of A(t) in the future happens much later. 

Models with the growing cosmological term energy density and the time-de- 
pendent G have testable predictions about the past evolution of the universe. 
The confrontation of these predictions with the host of cosmological observations 
available today gives constraints on the parameters of the model, above all on 
the parameter 77. These constraints determine the destiny of the universe and its 
structures and how far in the future the possible dramatic unbounding effects lie 
in front of us. 
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Figure 1: The evolution of the scale factor of the universe is shown for $7° = 0.7, 7 = 
and four different values of the parameter r\. The difference between various models is 
the most prominent in the distant past and the distant future. More negative values of 
rj lead to the faster growth of the scale factor in the distant future. 
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Figure 2: The evolution of the gravitational constant G is given for = 0.7, 7 = 
and four various values of the parameter 77. In the very early universe the change of G 
is very slow, while in the distant future G has a tendency towards 0. For more negative 
values of the parameter r], the quantity G changes faster. 
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Figure 3: The change in time of the CC is depicted for = 0.7, 7 = and four different 
values of the parameter rj. In the distant future the quantity A reaches a saturation value. 
The onset of saturation is sooner for the smaller values of the parameter rj. 
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Figure 4: The asymptotic value Aoo of the CC as a function of the parameter r\. Param- 
eters used in the calculation are il° A = 0.7, 7 = 0. For a displayed range of values of the 
parameter rj, the asymptotic value is of the same order of magnitude as the value of 
the CC today, A . 
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